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^ ! Abstract 

We investigate the asymptotic separation of the optical model potential for nucleon-nucleus 
scattering in momentum space, where the potential is split into a medium-independent term and 
another depending exclusively on the gradient of the density-dependent g matrix. This decomposi- 
tion confines the medium sensitivity of the nucleon-nucleus coupling to the surface of the nucleus. 
We examine this feature in the context of proton-nucleus scattering at beam energies between 30 
and 100 MeV and find that the pn coupling accounts for most of this sensitivity. Additionally, 
based on this general structure of the optical potential we are able to treat both, the medium de- 
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pendence of the effective interaction and the full mixed density as described by single-particle shell 
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models. The calculated scattering observables agree within 10% with those obtained by Arellano, 
^! ■ Brieva and Love in their momentum-space (/-folding approach. 
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I. INTRODUCTION 



Microscopic optical model potentials for nucleon-nucleus (NA) scattering are usually ex- 
pressed as the convolution of a two-body effective interaction with the target ground-state 
mixed density. Their realization becomes feasible with the use of Brueckner-Bethe-Goldstone 
g matrices and resorting to simplifying assumptions in their coordinate and/or momentum 
dependence. Thus, nuclear medium effects are disclosed by means of volume integrals of 
density-dependent interactions throughout the nucleus [l[ [2]. 

In this article we explore in more detail a recent finding suggesting that intrinsic nu- 
clear medium effects, namely those arising from the dependence of the g matrix on the 
density, are dominantly localized in the nuclear surface, i.e. regions where the gradient of 
g is strongest j^]. This result is a consequence of a close examination of the momentum- 
and coordinate-space structure of a two-body effective interaction spherically symmetric in 
its mean coordinate. It is demonstrated quite generally that two-body interactions can be 
expressed as a non-trivial sum of a medium-independent term and another which is function- 
ally -and exclusively- proportional to the gradient of the reduced in-medium interaction. 
As a result, the optical potential in momentum space becomes the sum of medium-free (free 
t matrix) and medium-dependent (g matrix) contributions, the latter depending exclusively 
on the variations of the effective interaction with respect to the mean coordinate. This 
feature yields an enhancement of the intrinsic medium effects in the nuclear surface and 
suppression in the saturated volume. We have investigated the manifestation of this selec- 
tivity in NA scattering at proton energies below 100 MeV, identifying major sensitivity in 
the pn couplings. 

Current trends in nuclear research and applications have resulted in the development and 
construction of novel research facilities around the world. Such is the case of radioactive ion 
beam accelerators in the US, Europe and Japan 0, [5, f], B, [^j], where intense rare isotope 
beams are produced and collided against selected targets. As the energy of these beams are 
projected to reach 500A MeV, their scattering from hydrogen targets would result equivalent 
to intermediate-energy proton collisions from an exotic nucleus. At energies below 60A MeV 
the collision would correspond to low-energy proton-nucleus scattering, a regime where the 
inclusion of medium effects in the effective interaction is known to be important. These 
unique facilities will expand our access to the neutron drip line from the region roughly 
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below carbon isotopes to nuclei as heavy as 52 S, the most neutron-rich nucleus. Just before 
the neutron drip line is reached, neutrons occupy weakly bound states spread well away 
from the bulk of the nucleus. Hence, from the NA scattering point of view, counting on an 
approach capable of tracking more selectively the various contributions to the optical model 
potential -particularly its surface structure- could prove useful for studying and interpreting 
data from rare isotope beam facilities. In this article we present and investigate a simplified 
form of the unabridged optical model potential discussed in Ref. j^], paying attention to 
the surface structure which emerges from the intrinsic medium effects as implied by the 
theory. Additionally, this form allows us to treat explicitly the off-shell mixed density, a 
long standing limitation of the microscopic in-medium folding approach of Arellano, Brieva 
and Love (ABL) [9j, where the Slater approximation is used. 

This article is organized as follows. In Section II we outline the general framework, discuss 
the general structure of the optical model potential in the single scattering approximation, 
introduce the '£g-folding' approach and make contact with known approximations. Addi- 
tionally, we examine more closely the various contributions to the optical potential, their 
energy and density dependence. In Section III we present and discuss results from selected 
applications of proton elastic scattering at energies between 30 and 100 MeV. In Section IV 
we summarize this work and present its main conclusions. 



II. THEORETICAL FRAMEWORK 



From a broad perspective diverse formal expressions of the optical model potential for 
NA scattering can be found in the literature 
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12L [13|, [14(. Although they may differ in 



the way contact is made with the bare NN potential, they all take the form of a ground-state 
expectation value of a generalized two-body interaction. Thus, a general representation of 
the optical model potential for collisions of nucleons with kinetic energy E off a composite 
target is given by 



U(k', k;E) = / dp' dp (k'p' | f(E) \kp) A p(p',p) , 



(1) 



where the subscript A indicates antisymmetrization. In general, T contains information 
about the discrete spectrum of the many-body system. In this expression p(p', p) represents 
the one-body mixed density corresponding to the ground-state of the target. A compre- 



3 



hensive evaluation of the optical potential considering the full f matrix would require the 
solution of the (A+ l)-body system, a formidable challenge. This difficulty is circumvented 
by decoupling the two-body effective interaction from the ground-state structure, a suitable 
strategy at intermediate and high energies when the discrete spectrum of the many-body 
Green's function is away from the projectile energy in the continuum. This allows the 
use of single-particle models to describe the target ground state and the Brueckner-Bethe- 
Goldstone reaction matrix to represent the effective interaction. 

A. Structure of the optical potential 

As expressed in Eq. ([TJ, a central element in the evaluation of the optical potential 
is the representation of the two-body effective interaction. Quite generally, regardless of 
the physics content or particular structure conceived for the NN interaction, the two-body 
operator T in coordinate space requires the specification of four vectors. This leads to matrix 
elements of the form (r's' | T \ r s), where r and s denote the 'prior' coordinates of the 
projectile and target nucleon, respectively. The primed vectors refer to 'post' coordinates. 
With these definitions, the so called local coordinate Z (hereafter referred as the mean 
coordinate) becomes 

Z = l(r' + r + s' + s) , 

corresponding to the simple average of the prior and post coordinates of the interacting pair. 

As demonstrated in Ref. the momentum-space representation of the T matrix can be 
cast in terms of a reduced interaction, gz, in the form 

(feV \f\kp) = J e iZ <^ g z (K li; k', k) . (2) 

In Eq. (T5]) Q = p' — p, the recoil of the target nucleon; q = k — k', the momentum transfer 
of the projectile; K = (fe + k')/2, the mean momentum of the projectile; P = (p' + p)/2, 
the mean momentum of the struck nucleon; and 

k' = \{k' - p') = \{K - P - | (q + Q)] , (3a) 
K=i(fc-p) = \[K-P + \{q + Q)\ , (3b) 

the post and prior relative momenta, respectively. Furthermore, 

K ll = K + P = l(k + k' + p + p') , (4) 
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interpreted as a longitudinal momentum of the interacting nucleons 3j. With these defini- 
tions the integrals on (p,p f ) in Eq. are accounted for by (P, Q), with dp' dp = dQdP. 

What is interesting about the above representation for T is that it prescribes naturally 
the way the medium dependence of the two-body interaction is mapped through the mean 
coordinate Z in the reduced interaction. To model this dependence we have resorted to 
infinite nuclear matter, a reasonable starting point to incorporate leading-order correlations 
in the nuclear medium. In this approach, to each coordinate Z we associate its nuclear 
isoscalar density p(Z) = [p n (Z)+p p (Z)]/2, therefore a symmetric nuclear matter Brueckner- 
Bethe-Goldstone reaction matrix gz {g matrix) satisfying 

^ =d + d n n Q l r m ■ (5) 

i I + iv — hi — h 2 

Here v is the bare NN potential, hi and h 2 the quasi-particle energies at density p, and Q 
the Pauli blocking operator to suppress occupied intermediate states. The corresponding 
Fermi momentum is given by 

k F = (3vr 2 p) 1/3 • (6) 

In a finite system, namely a system with confined matter distribution, we demand that 
lim^oo p(Z) = 0, so that 

lim g z (Q) = t(Q) , (7) 

Z— >oo 

the free-space t matrix. 

In the context of a spherically symmetric matter distribution, the Z integral in Eq. ([2]) can 
be split in such a way that its asymptotic structure becomes isolated from the Z-dependent 
term [3]. Accordingly 

(k'p' \f\kp) = 5(Q-q) t(K f k', ^-^J™ z3dz *i(Z\Q~Q\)^§, (8) 

where the momentum dependence of dgz/dZ on K\\, k,' and k, is implicit. Here $i(i) = 
ji(t)/t, with ji the spherical Bessel function of order 1. This profile function favors the 
recoil of the struck nucleon around Q ~ q, namely k + p w k' + p'. We observe that 
total momentum conservation can only be possible when the system exhibits translational 
invariance, as expressed when dgz/dZ = 0. 

Upon substitution of T from Eq. ([S]) into Eq. (pQ) for the optical potential we obtain 

U = U + U 1 , (9) 



with 

U {k',k;E) = J dP p(q;P)t(KfK',K) ; (10a) 
U x (k',k',E) = dQdPp(Q;P) x 

J™ Z'dZ <S> X {Z\Q - q\) (-^fj . (10b) 

The first term, Uq, depends exclusively on the medium-free reduced matrix, whereas the 
second depends on the gradient of the g matrix. In these expressions p denotes the full 
mixed-density, which in terms of occupied single-particle states <f> a is given by 

p(Q; P) = ^p Q (Q; P) = J2^(P') MP) ■ 

a a 

Eq. (Q for U represents the most general expression to be given to the optical model 
potential when the two-body effective interaction exhibits spherical symmetry in the mean 
coordinate Z. It summarizes the medium dependence of a general two-body effective inter- 
action, accounting for all phase-space configurations allowed by the one-body mixed density. 
The interaction is evaluated off-shell, with no assumptions regarding its local/nonlocal struc- 
ture. Furthermore, it involves a sevenfold integral, sixfold in momentum space and an ad- 
ditional integration in coordinate space. Thus, its evaluation constitutes a very challenging 
task even for nowadays computational capabilities. In this work we circumvent this difficulty 
by introducing a simplifying assumption, within the momentum- conserving approximation, 
to be explained in the following subsection. A thorough assessment of this assumption and 
its implications in actual scattering processes may require the evaluation of the sevenfold 
optical potential itself. 



B. Limit cases and the ^-folding 

The general form of the optical potential expressed above leads naturally to the free t 
matrix and ABL folding approaches. For instance, if the effective interaction is taken as 
the (transitionally invariant) free t matrix, then U\ vanishes and U becomes Uq, the free 
t-matrix full-folding optical model potential applied to intermediate-energy NA scattering 



15l . Il6| . In this case the medium effects do not come from the effective interaction but from 
the Fermi motion of the struck nucleons, as allowed by the spread of the one-body mixed 
density. 



The evaluation of the optical potential for a spherically symmetric system involves a 
7-fold integral. In order to simplify this, we neglect the dependence of dgz/dZ on Q by 
setting Q — > q in the interaction. This change is designated as gz — > g*§\ so that 

^«^Ef z3rfz / dp n ^ p > z ) (-if-) ' ( n ) 

with 

fl a (qf,P;Z) = J dQ p a (Q;P) <S> X {Z\Q - q\) . 

Noting that Q — > q expresses momentum conservation of the interacting pair in the g matrix 
(k+p = k'+p'), we find appropriate to refer to this as momentum-conserving approximation 
(MCA). The appealing feature of this result is that it enables a detailed treatment of the 
full-mixed density as obtained from single-particle shell models while accounting for the 
medium dependence in the g matrix. The following discussions are mainly focused on this 
structure of the optical potential, referred hereafter as { 8g- folding'. 

As demonstrated in Ref. [3], the use of the Slater approximation within the MCA leads 
to the ABL potential 

U ABL = 4vr / Z 2 dZ Jo (qZ) p(Z) (gf) , (12) 
Jo 

where (g^ ) denotes the Fermi-motion integral 

(gP) = J dP S F (P; Z) g z [K h \{K-P-q),\{K-P+ q)} , 

with 

s f (p ] z) = t J— e[k(z)-p]. 

This step-function sets bounds for the off-shell sampling of the g matrix at a distance Z 
from the center of the nucleus, | P \< k(Z), with k(Z) = [3ir 2 p(Z)] 1 ^ 3 . The above result for 
Uabl can also be obtained by replacing directly the two-body effective interaction [c.f. Eq. 
(T5])] into Eq. ([I]) for U, applying the MCA and representing the mixed density by its Slater 
form. 

All the above forms of the optical potential are nonlocal, as a consequence of the momen- 
tum structure of the g matrix -solution of the Brueckner-Bethe-Goldstone integral equation- 
expressed in terms of the relative momenta n 1 and k [c.f. Eqs. ( I3a|3bj) ]. The antisymmetriza- 
tion of the interaction accounts for additional nonlocalities. These features have not been 
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duly explained in previous works, leaving room for misconceptions. So it may be worth 
to sketch them here for clarity. To make the illustration simple let us consider the rank-0 
(scalar) antisymmetrized reduced g matrix for total spin S and isospin T, 

(«' I 9 ST I *U = 0*V, k) - R S+ VV, -k) . 
Expanding in partial waves 

L=0 

and using the property Pl(—u) = (—) l Pl(u), we can arrange the antisymmetrized g in a 
single sum, 

(K> | g ST | k) a = £>£V, k) [1 - P L {k' ■ k) . 

L=0 

Therefore 

(*> | g ST | k>^ = 2 ^ «) P l (k' ■ k) , (13) 

Allowed 

where the summation considers only those NN states allowed by the Pauli exclusion prin- 
ciple and the off-shell matrix elements g£ T (ft', k), direct solutions to the Brueckner-Bethe- 
Goldstone equation for the corresponding partial wave. In this fashion we naturally account 
for the knock-out exchange term. 

In the case of local effective interactions , the off-shell matrix elements q (k', k) 

are obtained via the Fourier transform v of the local function v ST (r), hence 

/ T « k) = v st (k> - K ) . 

Therefore, the antisymmetrized matrix element reads 

(«' | g ST | k)a = v st (k' - k) - (-) s+t v st (k' + k) , 

a well known result. Here the knock-out exchange term makes the antisymmetrized inter- 
action non local. If one uses multipole expansions to these Fourier transforms, then the 
antisymmetrized {g ST )A takes the same form as that expressed by Eq. ( |T3l) . In this case, 
however, <?£ T (/c', k) = r 2 j L (n'r)v ST (r) j L (nr) dr. 
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C. The medium-dependent term 



We examine more closely the structure of Ui, particularly the shape of its integrands. 
Since the dependence of g matrix elements on Z is set via the isoscalar density p, with 
P = kp/37r 2 , then we can write 



In Fig. ([I]) we plot the radial dependence of the density p(Z) (upper frame), its cor- 



—k'(Z) (lower frame), for 16 (solid curves) and 90 Zr (dashed curves), respectively. These 
figures exhibit clear peaks of —k'(Z) near 3 fm and 6 fm, corresponding in both cases to k ~ 
0.6 fm -1 , i.e. the surface of the nucleus. We estimate in ~3 fm the width of both peaks, 
limiting the region where the main contributions to U\ should come from. The strength 
of these contributions are dictated by the derivative Sg = dg/dkp, which depends on the 
energy E of the projectile. 

In Fig. (j2J) we show the partial derivative of the on-shell g amplitude with respect to 
the Fermi momentum, symbolized with 5g. The real and imaginary components are shown 
in the upper and lower frames, respectively. The left frames correspond to the pp channel, 
whereas the right frames correspond to the pn channel. The curves represent different 
projectile energies, starting at E =30 MeV (solid curves) up to 100 MeV in steps of 10 MeV 
(dashed curves). To facilitate their comparison, the same scale is used in all graphs. By 
forward (on shell) we mean k' = k, with E = k 2 /2m, the nonrelativistic nucleon energy. 
Hence, 5g = dg(k; |fc, ^k)/dkp. The striking feature of this figure is the asymmetry of 6g, 
significantly more pronounced in the pn than in the pp channel, suggesting more sensitivity 
to neutron densities when protons are used as projectiles. Looking at the real part of the pn 
coupling, the attraction is more pronounced in the region 0.2 fm -1 < kp < 0.6 fm -1 , i.e. the 
nuclear surface, a feature which diminishes with increasing energy. Regarding the imaginary 
contribution, the nuclear surface contributes with more absorption, whereas in the nuclear 
interior (kE > 0.6 fm -1 ) it is weakened. It is important in this analysis to keep in mind 
that 5g contributes to U\. Instead, the leading-order contribution to the optical potential 




with 



k\Z) 



k(Z) dhxp 



3 dZ 



responding local Fermi momentum k(Z) (middle frame) and the negative radial derivative 
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stems from Uq, which depends directly on the t matrix. To keep this observations in better 
perspective, in Fig. ([3]) we plot the forward (on shell) t matrix as function of the nucleon 
energy E. Here the solid and dashed curves correspond to the real and imaginary amplitude, 
respectively. In this figure the right-hand-side axis scales to (27r) 3 t, to facilitate comparison 
with other conventional normalizations. We notice here that the absorptive component of 
the pn coupling exhibits a stronger energy dependence, becoming dominant as the energy 
decreases from < 80 MeV. Instead, the real components are relatively constant throughout 
the energy range considered. 

In order to trace the sources of the contributions to U\ and also estimate their importance 
relative to U , we find useful to introduce the density function u a defined by 



to be evaluated on-shell at q = for the single-particle shell a. Its radial integral accounts 
for the partial contributioin U a , 



with Ui = Yl a Ua- In Fig- © we plot the real (upper frames) and imaginary (lower frames) 
components of u a (Z) for w O(p,p) at 30 MeV. The curves correspond to contributions from 
the lpz/2 (dotted curves), lpi/2 (dashed curves) and lsi/ 2 (long-dashed curves) shells, while 
the solid curves represent the sum ^2 a Ua' nS) - The (p) and (n) labels symbolize contribu- 
tions of the form {p p Sg pp ), arising from proton densities, and (p n 5g pn ) due to neutrons, 
respectively. Notice that the scale of Im u a doubles that of Re u a . 

This figure evidences quite neatly surface-peaked structure stemming from 5g, confined in 
the region 3-5.5 fm, with clear dominance of neutron over proton distributions. Considering 
Uq the leading-order contribution, Re Ua enhances the attraction to the projectile (proton). 
This can be readily estimated considering its width ~1.5 fm and depth ~15 MeV fm 2 . 
Hence, the area between the curve and the Z axis yields Re U[ n) ~ -23 MeV fm 3 . This 
is to be compared to ~ 8xRe t pn ~ -10 MeV fm 3 , as extracted from Fig. ([3]). In 
contrast, neutron density contributions to Im U\ becomes considerably weaker due to its 
near-canceling up-and-down structure observed in the lower-right frame, while Im XJ^ ~ 
8xlm t pn ~ -24 MeV fm 3 . The extent to which these features become important in collision 
processes needs to be assessed by examining scattering observables. In any case, the pocket 




(14) 



a 
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shape of Re Ua near the surface indicates a preference to couple the projectile (proton) 
with the ulpi/2 and z/lp 3 / 2 shells, favoring (p, d) pickup reactions. This feature is consistent 
with recent findings on pickup effects in p+ 10 Be elastic scattering near 40A MeV frol ]. 



III. APPLICATIONS 

We investigate proton elastic scattering from 16 and 90 Zr, two relatively well known 
doubly closed-shell nuclei. In each case we consider three forms of the optical potential. First, 
the Sg- folding approach [c.f. Eq. (1111) ]. providing arguably the most complete momentum- 
space description of the optical potential to date. Here, single-particle wavefunctions are 
used to represent the one-body full mixed density while a thorough account of the medium 
dependence of the antisymmetrized off shell g matrix is given. Second, the ABL folding 
approach, corresponding to a simplified representation of the mixed density in terms of its 
Slater approximation. This approach has been extensively discussed in Refs. jsl, Q] . Lastly, 



the free t-matrix full-folding optica 



used as in the early calculations 
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potentials (t- folding) , where the full mixed density is 



id 



The calculations presented here are based on the Paris NN potential [21[ . We have inves- 
tigated other NN potentials and found no significant differences with the results reported 
here. The corresponding g matrix was calculated off shell ( J < 7) at 30 values of the Fermi 
momentum, ranging from up to 1.6 fm -1 . This thin mesh is no longer necessary after 
various tests of convergence were performed; the use of 16 Fermi momenta yields equally 
reliable results. The needed selfconsistent nuclear-matter fields were computed prior to all 
runs. 

To evaluate U\ given by Eq. (TTTj) we carry out the P and Q integrals using Gauss- 
Laguerre quadrature at 25 radial mesh points. The Z integration is performed using a 
uniform mesh with steps of 0.1 fm. As stated earlier, the calculated optical model potentials 
reported here are nonlocal operators, treated as such throughout. The scattering observables 
are obtained solving the Schrodinger equation with the nonlocal coupling in the presence of 
the Coulomb term. See Ref. 22 1 for more details. 
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A. p+ Zr scattering 

In Fig. (jSJ) we present the measured and calculated differential cross section, as a function 
of the center-of-mass scattering angle 6 c . m ., for 90 Zi(p,p) at 30.4 and 40 MeV. The data are 



from Refs. [23| and [24J, respectively. The solid curves correspond to Sg- folding, the dashed 
curves to the ABL approach and the dotted curves to the t-folding. The one-body mixed 
density is constructed using single-particle wavef unctions based on Hartree-Fock calculations 



27] . 



The calculated cross sections based on the 5g-folding follow reasonably well the diffractive 
pattern exhibited by the measurements. The maxima are in phase with the data, although 
the diffractive minima tend to be more pronounced. Additionally, we find that the ABL 
approach (dashed curves) follows very closely the Sg- folding (solid curves). So, within the 
MCA, the ABL folding approach represents reasonably well the <5g-folding. The difference 
lies in the computational time needed for their evaluation, being the <5g-folding more time 
demanding (by nearly a factor 150) than its ABL counterpart. 

Considering the results based on the free t matrix (dotted curves), they clearly lack the 
structure exhibited by the data. In particular, at the two energies considered here, the 
first diffractive minima occur at greater angles than those shown by the data, suggesting a 
smaller size nucleus. Thus, the medium effects accounted for by Ui, and located mainly in 
the nuclear surface, do account for some of the hadronic size of the nucleus. 

As the energy increases it is expected that the medium effects become less relevant in the 
scattering process. This feature is clearly observed in Fig. (jSJ), where we plot the measured 



and calcu 
Refs 



ated differential cross section for 90 Zr(p,p) at 80 and 100 MeV. The data are from 



251 ] and 26], respectively. The curves follow the same convention as in Fig, (j3j). 
The agreement of the 5g-folding with the data is remarkable throughout the whole range of 
the measurements. Additionally, we verify that the difference between the t-folding results 
and those from 5g- or ABL-approach have diminished considerably relative the previous 
applications. 
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B. p+ O scattering 

In Fig. ([7]) we show the measured and calculated differential cross section (upper frames) 
and analyzing power (lower frames), as a functions of the center-of-mass scattering angle, 



for w O(p,p) at 30.4 and 49.48 MeV. The cross section data are from Refs. [28j and 29], 
respectively. Here again we observe that the 5g-folding and ABL approach follow very close 
each other. Also, the free t matrix results agree poorly with the data, as expected. This 
lighter target evidences a deviation of the calculated cross sections with the data, particularly 
the depth of the diffractive minima in the cross section. Indeed, both 5g- and ABL-folding 
approaches fail to account for the shallow first minima near 40°. In turn, both yield non 
existing minima near 90° and 80°, respectively. 

We have performed various tests of sensitivity to assess the consistency of the results 
presented here. For instance, using harmonic-oscillator wavefunctions for 16 -with the 
same root-mean-squared radii- we obtain practically the same results for the scattering 
observables. To keep this work focused on the structure of the 5g-folding, we have not 
explored the sensitivity of the calculated scattering to alternative representations of the 
target ground state, leaving such study for future works. 

In Fig. (jSJ) we examine how the differential cross section is affected upon changes on 
U\. The case is 16 0(p,p) at 30.4 MeV. Here the solid curve represents results from the 
(^-folding, Uq + Ui, while the dotted curve is based on U alone, namely the t-folding. 
What is interesting to note is the effect of suppressing selectively the proton density (dash- 
dotted curve) and neutron density contributions to U\ (dashed curve). We note that the 
role of is considerably more significant than that of . Indeed, by neglecting the 
cross section ends up being very similar to the one obtained with the t-folding, in contrast 
with the moderate change on the 5^-folding result when Z7j is suppressed. These results 
are consistent with our analysis of u a (Z) discussed in the previous section, confirming the 
importance of neutron distributions in the optical model potential at these low energies. 

IV. SUMMARY AND CONCLUSIONS 

We have investigated the structure of the optical model potential as inferred from its 
general form, once the MCA is applied to the vector structure of the NN couplings. The 
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resulting (^/-folding) potential becomes expressed as the sum of two components, Uq + 
Ui, where Uq corresponds to the free t matrix full-folding potential and U\ folds the full 
mixed density with the gradient of the medium- dependent effective interaction. This feature 
implies that the intrinsic medium effects are localized mainly in the nuclear surface. When 
comparing the relative strength of these contributions, we find that the pn coupling is 
considerably stronger than its pp counterpart, a feature that fades out as the energy of 
the projectile is increased. This asymmetry leads to stronger medium-sensitivity of proton 
scattering to neutron matter distributions of the nucleus. 

As a by product of this study, with the introduction of the <5g-folding we have been able to 
provide a thorough account of the full mixed density in the evaluation of momentum-space 
optical potentials within the g matrix. With this we mend a long standing limitation of the 
ABL folding approach to NA scattering, where the mixed density has been approximated 
by its Slater form. When comparing the differential cross section, the <5g-folding and ABL 
approach are close to one another within 10% in the diffractive maxima. 

We have assessed the predicting power of the <5g-folding approach to proton elastic scat- 
tering from 16 and 90 Zr at energies between 30 and 100 MeV. In the case of 90 Zr(p,p) 
we are able to provide reasonably good descriptions of the data. For 16 0(p,p), in turn, 
the differential cross section is underestimated significantly. At this point we are not clear 
whether these shortcomings stem from missing contributions implied when the dependence 
of dg/dZ on Q is neglected, the existence of exotic neutron structures in the surface, or from 
higher-order effects in the in-medium effective interaction. One has to keep in mind that the 
low-energy interaction of the projectile with target nucleons becomes more sensitive to the 
shell structure of the nucleus, in addition to the presence of collective excitations or other 
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30]. 



reaction channels. Investigations along these lines have recently been reported 

The present work constitutes a step forward toward a comprehensive momentum-space 
description of the optical model potential for NA scattering, in the form of a unified descrip- 
tion, for elastic and inelastic processes from few tens of MeV up to GeV energies. The only 
microscopic inputs to achieve this goal are the bare NN potential and the target ground-state 
mixed density, although high-energy applications may also require NN phase-shift analyzes 
accounting for loss of flux above pion-production threshold |20j. The introduction of the 
5g-folding optical potential has allowed us to visualize very simply the interplay among dif- 
ferent elements in the interaction of a single nucleon with finite nuclei, particularly the role 
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of the in medium interaction in the nuclear surface. With the application reported here we 
have been able to set narrower margins of uncertainty in the evaluation of the first-order 
optical model potential, an important consideration for high-precision analyzes of upcoming 
scattering data involving unstable nuclei. 
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FIG. 1: Radial dependence of the density (upper frame), local Fermi momentum (middle frame) 
and its negative gradient (lower frame) for 16 O and 90 Zr. 
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FIG. 2: Real (upper frames) and imaginary (lower frames) components of the forward 5g = dg/dkp 
amplitudes, in the pp (left frames) and pn (right frames) channels, as functions of the Fermi 
momentum. See the text for reference to the curves. 
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FIG. 3: Forward on-shell free t matrix as function of the nucleon energy. The solid and dashed 
curves correspond to the real and imaginary amplitudes, respectively. 
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FIG. 4: The radial behavior of u a for the pp (left frames) and pn couplings (right frame). Observe 
that J2 a I uS ,n \z) dZ = Jj[ p ' n \ See text for reference to the curves. 
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FIG. 6: The measured and calculated differential cross section for 90 Zr(»,w) at 80 and 100 MeV. 

I n 

The data are from Refs. [25[ and [26], respectively. See text for reference to the curves. 
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FIG. 7: Measured and calculated differential cross section for 1& 0(p,p) at 30.4 and 49.48 MeV 
(upper frames), and calculated analyzing power (lower frames). The data are from Refs. y and 



29J, respectively. See text for reference to the curves. 
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FIG. 8: The calculated differential cross section based on the 5q- folding (solid curve) and t-folding 



(dotted curve) for 16 0(p,p) scattering at 30.4 MeV. The data are from 
to the dashed curves. 



281 ] . See text for reference 
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